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Abstract 


This study identified the fluid types and boundaries present within selected reservoirs in 
Langbodo field, using petrophysical parameters based on estimated rock properties such as 
porosity, permeability, irreducible water saturation, hydrocarbon saturation and bulk water 
volume. This was with a view to correcting the salient reservoirs heterogeneities anomalies 
error inherent in building of an ideal realistic reservoir models. The quality of the data obtained 
were checked and despiked to eliminate null values. Petrel version 2009 and OpendTect 4.6.0 
Exploration and production softwares were used for the quality interpretations of data, such as 
lithology identification, delineation of potential reservoirs and determination of fluids and 
fluids contacts. Estimation of quantitative petrophysical parameters were done by inputting the 
data into Microsoft excel 2015 version softwares and adopting appropriate mathematical 
relations, such as the Tixier, Timur and the Coates and Dumanoir models for the permeability 
(K). Realistic estimation of the permeability was done by comparing the average of the Tixier, 
Timur and the Coates and Dumanoir models with each of the models. The composite model 
obtained, mirrors the behavior of the Timur’s permeability which is higher than that of the 
Tixier and the Coates and Dumanoir. Integration of the Achie’s equation and neutron — density 
crossplot confirmed the presence of substantial hydrocarbon in the reservoirs, although 
producibility indicators revealed that the reservoirs may not be producible without enhanced 
oil recovery method(s). This study established that the composite model is a better 
representation of K in the study area because it agrees with the Timur’s estimation model. 
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I. INTRODUCTION 


he geometric progression of the world population has 

called for higher demand for energy, of which 
hydrocarbons constitute a dominant percentage, especially 
among the non-renewable energy sources. Multinational 
hydrocarbon exploration companies may experience poor 
reservoir performance within few years of production due to 
inadequate reservoir properties description [1]. The success of 
any hydrocarbon exploration program depends on the building 
of reliable reservoir models [1]. Some of the most important 
parameters needed to characterize reservoir quality are 
permeability, porosity and shale volume. Their accurate 
prediction is the basis upon which one can actually identify if 
the reservoir is producible. The description of reservoir 
characteristics and fluid flow performance can be anchored on 
the permeability and this plays a very important role in 
designing exploration and development plan [2]. An accurate 
estimation of the permeability enhances oil and gas field 
development and it is the basis for building geophysical 
models, accurately predicting oil and gas reserves and taking 
reasonable development plan [3]. The importance of accurate 
estimation of the permeability cannot be overemphasized as 
other petrophysical parameters are dependent directly or 
indirectly on it. Therefore, accurate estimation of the 
permeability in identified reservoirs is undoubtedly important 
to enhance validity of other dependent variables. 
Unfortunately, diverse methods, with varying strengths and 
weaknesses have been proposed by various researchers. The 
Tixier’s model of 1949, the Timur’s model of 1968 and the 
Coates and Dumanoir’s model of 1981 [4-6], were applied by 
several researchers [7-10], with little or no cognizance to the 
weaknesses inherent in such approaches. These models are 
based on correlation between permeability and other 
geophysical parameters such as porosity and irreducible water 
saturation. The estimation of permeability from the 
aforementioned models is often accompanied with varying 
strengths and weaknesses; which the ignorant and 
unsuspecting explorationist may not watch out for; depending 
on the peculiarity of the geology of the area under 
consideration. Therefore, this research work is aimed at 
enhancing the computation of the permeability by integrating 
the strengths of the diverse methods for estimating the 
parameter, using suites of well logs from “Langbodo field” 
Niger Delta. This research finding will improve reservoir 
quality assessment, and assist other producibility indicators to 
rank the reservoirs for further developmental decisions. 


Il. STUDY AREA 
The study area is located in the “Langbodo field” onshore 
Niger Delta. It is located within the Shell Petroleum 
Development Company acreage. The field is bounded by 
latitudes 4°46’N and 5°57’N and longitudes 5°37'E 
and 5° 64’ E. Fig. 1 and 2 are maps showing the study area 
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and well location points respectively. The subsurface geology 
of the study area reveals the Niger Delta Basin. The wells 
drilled in the study area enabled the acquisition of wireline 
logs, utilized for this study. 


Fig. 1 An Outlay of the “Langbodo Field” showing Well 
Locations 


The geology of the Niger Delta revealed an extensional rift 
basin located in the Niger Delta and the Gulf of Guinea on 
continental margin close to the western coast of Nigeria. It has 
an established access to Equatorial Guinea, Cameroon and Sao 
Tome and Principe. The complexity of the basin is shown in 
the content of high productive hydrocarbon system. The Niger 
Delta basin is one the largest sub-aerial basins in Africa which 
is composed of several geologic formations that indicate how 
the basin could have initially formed and the large scale 
tectonic of the area. Research shows that some other basins 
formed from similar geologic processes exist around the Niger 
Delta. Its formation can be traced to a failed rift junction 
formed during the separation of the South American plate and 
the African plate, as the South Atlantic started to open [10]. 
The origin of rifting in this basin started in the late Jurassic 
and terminated in the mid Cretaceous. Continuation of rifting 
led to the formation of several faults many of which are thrust 
faults. Also at this time, the late Cretaceous harbors some 
deposition majorly composed of syn-rift sand and shale. This 
shows that the shoreline regressed during this time. 
Concurrently, the basin had been witnessing extension leading 
to high angle normal faults and fault block rotation. At the 
origin of the Paleocene, there was a noticeable shoreline 
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transgression. During the Paleocene, the Akata Formation was 
deposited. At the time of the Eocene, the Agbada Formation 
followed the underlying Akata shale [12]. This formation 
loading caused the underlying shale Akata Formation to be 


compressed and squeezed into shale diapirs. The Oligocene, 
the Benin Formation was thereafter deposited. It is the 
shallowest part of the sequence with age of formation varying 
from earlier to recent [12]. 


Fig. 2 Location Map of the Petroliferous Niger Delta, Showing the Important Frame and Tertiary Delta Growth [11]. 


Ill. MATERIALS AND METHODS 


If Suites of wireline logs that include the gamma ray, 
density, neutron, sonic, resistivity and caliper logs from three 
exploration wells were utilized for this study. These data were 
checked and despiked to eliminate null values. Data were 
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thereafter loaded into petrel version 2009 and OpendTect 
4.6.0® exploration and production software to enable 
qualitative interpretations such as lithology identification, 
delineation of potential reservoirs, determination of reservoirs 
fluids and fluid contacts. Estimation of quantitative 
petrophysical parameters was done by inputting the data into 
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Microsoft excel® software and adopting appropriate 
mathematical/ statistical valuations to evaluate the delineated 
reservoirs. 

Enhanced computation of the permeability (K) of 
reservoirs involve a detailed qualitative and quantitative 
evaluation of other petrophysical parameters such as the 
lithology, porosity, water and hydrocarbon saturation 
among others. The analytic procedure is aimed at a better 
estimation of permeability (K). This will as well assist in 
determining quality of the reservoirs. Results obtained 
from at least three methods of estimating the reservoir 
permeability were used to arrive at a better K of the 
reservoirs. A systemic analysis suitable for this case was 
adopted for the estimation of the permeability of reservoirs 
of interest penetrated by the “Langbodo” wells. The 
Schlumberger Petrel Software was used to enhance 
interpretation. Microsoft excel was used to estimate K 
from various equations; as well as comparing obtained 
results. Table I below summarizes the data set used for this 
research. The Wells A, B and C penetrated a total 
of 11500 ft, 11620 ft and 12,035 ft respectively. 


Table I Available logs from the three wells used for the 
research 
Wells Gamma Density Neutron Resistivity Spontaneous Sonic RT Caliper 


ray potential 
A A A A A A A N A 
B A A A A N A A A 
C A A A A A A N A 
Legend 
A Available 


N Not Available 


IV. DETERMINATION OF RESERVOIRS 
PERMEABILITY 


Permeability is the rock property to transmit fluids. It is 
determined by the size of the connecting passages (pores 
throats or capillaries) between pores. Permeability is a key 
parameter associated with the characterization of any 
hydrocarbon reservoir. It is measured in Darcies or 
millidarcies. The Tixier (1949), Timur (1968) and the Coates 
and Dumanoir (1981) equations were used to derive the 
permeability of each reservoir that was identified on the 
Langbodo field. 


A. Permeability Estimation by Tixier Method 


Tixier, 1949 proposed that permeability (K) is a direct 
function of reservoir cubic porosity and an indirect function of 
the reservoir irreducible water saturation. Tixier’s equation is 
mathematically documented as follows: 

p? 
Krixier = 250 


Swirr 


(1) 


Where Krixier = Tixier Permeability, p= Porosity Swirr = 


Irreducible water saturation. 
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B. Permeability Estimation by Timur Method 

The second permeability equation used for this comparative 
study was proposed by Timur (1968). This equation is based 
on similar algorithm as that of Tixier (1949); although some 
variations exist between the duo. The earlier estimated 
irreducible water saturation was utilized to compute Timur’s 
permeability across 0.5 feet thickness of the bore wells; 
average values across each of the reservoirs were then taken. 


(2) 


44 
Krimar = 2505 


wirr 


Where Krimar = Timur Permeability, p= Porosity, Swirr = 
Irreducible water saturation. 


C. Permeability Estimation by Coates and Dumanoir 
Method 

Thirteen years after the emergence of Timur’s relation, Coates 
and Dumanoir proposed another relation that expressed the 
reservoir permeability as a direct function of the square 
porosity and a fraction of the irreducible water saturation. 


(3) 


Where Kcp = Coates and Dumanoir Permeability, = 
Porosity, Swirr = Irreducible water saturation. 


Kop = 100 9” (1-Swirr) 


Swirr 


V. DETERMINATION OF POROSITY 


Porosity is the ratio of the total volume of pore space to the 
entire volume of the formation. It relates the amount of 
internal space in a given volume of rock to the total volume of 
rock in the reservoirs. The amount of internal space or voids 
in a given volume of rock determines the amount of fluid it 
will hold. 


volume of pore space 


Total Porosity(%) = x 100% (4) 


total volume of formation 
Porosity may be classified as primary or inter-granular and 
secondary porosity. Primary porosity is the porosity that has 
existed in the formation since the time they were deposited. 
Secondary porosity is as a result of the action of tectonic 
forces or formation water. Porosity may be classified as 
effective and absolute porosity. Effective porosity is the 
porosity available to free fluid excluding unconnected 
porosity occupied by water bound and disseminated shale. 
Absolute porosity is the total porosity regardless of whether or 
not the individual voids are connected. Porosity was estimated 
from the density logs in this study. 


A. Density Log Derived Porosity 

Density log is a porosity log that measures electron density of 
a formation. It can assist the geologist to identify minerals, 
detect gas-bearing zones, determine hydrocarbon density and 
evaluated shaly sand reservoirs and complex lithologies. 


(5) 


Pma- Pb 


Porosity(g) = rar 
ma Pf 
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Where Pma = matrix density, pp =Formation’s bulk density 
and py = formation’s fluid density. 


Effective porosity was calculated for evaluating interval using 
the equation shown below: 


(@e = ema Pb _ Ven 


Pma- Pb 
Pma- Pf { } 


Pma- Pf 


(6) 


Where Pma= matrix density (usually 2.66 g/ce sandstone), 
Pp= Formation’s bulk density (obtained from density log at 
0.6ft. interval), p= formation’s fluid density (1.5 gm/cc for 
water and 0.8g/cc for hydrocarbon) and V,,= Density of 
adjacent shale body. 


B. Determination of Irreducible Water Saturation (S wirr) 
Irreducible water saturation is the water harbored in the pore 
spaces by capillary force. When a zone is at irreducible water 
saturation (Swirr), the water saturation in the uninvaded zone 
(Sy) will not migrate because it is held by pressure in the 
grains. For most reservoir rocks in the field, irreducible water 
saturation ranges from values less than 10 % to values more 
than 50 % [13]. The (Swirr) was estimated from the equation 
below: 


F 
Swirr = |= 
2000 


Where F = Formation factor and Swirr= Irreducible water 
saturation. 


(7) 


The formation factor was estimated from the Achie’s equation 
below: 


a 
F= om 


(8) 


Where ø = Porosity, a = Lithologic constant and M = 
Cementation exponent. 


VI. RESULTS AND DISCUSSION 


A. Quantitative Analysis of Reservoir Permeability (K) 


Recorded average permeability across the fifteen reservoirs 
mapped out are shown in Tables IHV. 


Table II Average Porosity and Permeability Values Obtained 
for Well A Reservoirs 


Porosity (p) Permeability (K) 
RSV Total Pe NPH F Ry Sw Krixier Krimur Kop Kaverage 
g (wr) I (Am) (vh) 


(viv) (PU) 
RSV1 0.316 0.289 0.348 16.593 1.490 0.505 244056.900 86717.010 3529.742 11434.500 
RSV2 0.277 0.252 0.339 17.412 1.567 0.795 16549.980 9091.426 979.819 8873.743 
RSV3 0.271 0.236 0.337 21.269 0.191 0.461 75419.440 24987.480 1183.889 33863.600 
RSV4 0.261 0.234 0.323 21.946 0.197 0.564 6088.246 4095.484 654.461 3612.720 
RSV5 0.269 0.237 0.328 21.081 0.189 0.842 2885.773 2955.228 646.456 2160.486 
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Table III Average Porosity and Permeability Values Obtained 
for Well B Reservoirs 


Porosity (p) 
RSV 


Permeability (K) 


Total 9, 
ọ (wv) I 
(v/v) (PU) 
0.256 0.246 0.366 
2 0.246 0.231 0.301 
0.234 0.222 0.288 21.013 0.168 0.271 
0.227 0.212 0.298 29.138 0.233 0.532 
5 0.219 0.204 0.277 1898.584 15.188 0.780 


NPH F Ry Sw 
(Am) (viv) 


Krixier Krimur Kep Kaverage 


18.895 
20.334 


0.151 
0.162 


0.166 
0.305 


57.378 
41.883 
29.210 
27.598 
23.398 


7075.634 71.716 
4292.459 56.940 
1388.265 45.814 
1395.567 42.829 
1065.497 37.746 


2401.576 
1463.761 
487.763 
488.664 
375.547 


Table IV Average Porosity and Permeability Values Obtained 
for Well C Reservoirs 


Porosity (p) 


RSV Total gy, 
g (vy) I 
(wv) (PU) 


Permeability (K) 
Krim ur Kep 


NPH F Re Sy 
(Qm) (viv) 


Krixier Kaverage 


RSV 
RSV2 
RSV3 
RSV4 
RSVS 


0.279 0.277 0.331 13.439 0.941 0.781 
0.275 0.267 0.320 17.453 1.221 0.980 
0.279 0.278 0.322 13.946 0.976 0.848 
0.274 0.271 0.320 14.941 1.311 0.987 
0.274 0.271 0.329 13.955 0.976 0.961 


8.259 
7.157 
8.269 
7.931 
7.891 


5810.766 9143.295 
5460.227 8620.817 
5757.388 9075.385 
5099.150 8085.388 
4954.188 7865.795 


4987.440 
4314.287 
4947.014 
4397.490 
4275.958 


Permeability obtained in well A by Tixier method ranged from 
244056.9 in Reservoir A to 2885.8 mD in Reservoir E. 
Similarly, Timur equation and Coates and Dumanoir equation 
yielded values that ranges from 86717.0 to 2955.2 mD and 
3529.742 to 640.5 mD respectively. Generally, permeability 
in well A is observed to decrease with depth. This may not be 
unconnected with close grain parking owing to overburden 
loads. Tixier method was observed to yield highest 
permeability across the reservoirs of well A. This trend 
amazingly change in well B, with Tixier recording values less 
than 100 mD, although these values were observed to be very 
close in range to the results obtained from Coates and 
Dumanoir equation. Based on the findings from the Langbodo 
field, the Tixier method was also observed to yield 
ridiculously low permeability values in interval penetrated by 
reservoirs (A to E) of “Langbodo” well C. The mean 
permeability values across the reservoirs are tabulated along 
the triple set of permeability equations Fig. 2, 3 and 4 are 
graphical representation of the triple set of equation and the 
mean value. The inconsistency noticed in the values obtained 
from Tixier method is believed to be due to certain effects on 
the reservoirs. This may be hydrocarbon effect or probably 
presence of heterogeneities that masked through permeability 
as presented by the method. [14], [15] and [16], have variously 
recorded permeability for some Niger Delta reservoirs: and the 
values obtained by the author are far more than that obtained 
by the Tixier method in wells A and B. Since mean values are 
true representations of dataset, the average permeability 
values are preferred to the other methods. It is however 
noticed that close similarities exist between the average 
permeability value and the Timur’s permeability (Fig. 3 — 5). 


©DOP_KASU Publishing 


PHYSICSAccess 


Damilare et al. 


Fig. 5 shows the comparison of the average permeability 
values to each of the values obtained from the triple set of 
equation across well A reservoirs. 


—=K Tixier 
—-—K Timur 
=+- CDK 

=r- Composite K 


Permeability ( 


1 2 3 4 5 1 2 3 4 


Reservoirs 
(a) (b) 


Fig. 3 (a and b) Comparison of Different Permeability 
Values Across Langbodo Well A 


The reservoir permeability obtained from RSV 1 of well A 
dropped suddenly in RSV 2. It increased a little in RSV 3 and 
later dropped in RSV 4 and 5. This result might be attributed 
to predominant laminated shale deposits in reservoirs 2, 4 and 
5 of well A. However, the Timur’s estimation model gave 
almost the same values with the Composite model which is a 
merger of the three other models. As earlier reported in this 
work, these discrepancies may be as a result certain 
heterogeneities such as shale/clay which masked through the 
permeability [16]. 


EEK Tixier 


E K Timur 
mmm coK 
[a Composte K 


Permeability (K) 


4 


Reservoirs 


Reservoirs 
ta) (b) 


Fig. 4 (a and b) Comparison of Different Permeability 
Values Across Langbodo Well B 


In well B, the Composite model was found to almost behave 
similar to the Timur’s estimation model. From this research 
we find higher estimated transmitivities and high reservoir 
permeability distribution within the models. The result from 
the Composite model shows a better representation of the 
permeability. This is similar to the work of [10], done by 
merging different models in the absence of core data from the 
study area to a nearby field core data in developing a reliable 
reservoir geophysical model. The differences in the behaviors 
of each of the models may be as a result of the permeability 
being masked by certain unnoticed heterogeneities [15]. 
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Permeability (K) 
Permeability (K) 


3 4 
Reservoirs 


Reservoirs 


Fig. 5 (a and b) Comparison of Different Permeability 
Values Across Langbodo Well C 


The three estimation models in well C for all the reservoirs 
behaved in almost similar ways. This may be as a result of 
similar deposits in each of the five reservoirs. It is evident 
from the Fig. 5 (a and b) above that the composite permeability 
estimation model which is obtained from the merger of the 
other three models employed in this work mirrors the Timur’s 
permeability model. Integration of the existing models to 
obtain a better reservoir permeability will give a better picture 
of the reservoirs in well C. This is similar to the work of [10]. 


B. Discussion of Other Petrophysical Parameters 


Some other petrophysical parameters were evaluated, part of 
the results obtained from these are shown in Tables 4.1, 4.2 
and 4.3. Other analyzed parameters such as the irreducible 
water saturation is tabulated in Table 4.4. Some of the 
important parameters evaluated are discussed below: 


1) Porosity 
Average values of porosity (Table II) obtained from the 
reservoirs of well A are, 0.3160, 0.2770, 0.2710, 0.2610 
and 0.2690 for RSV 1, RSV 2, RSV 3, RSV 4 and RSV 5 
respectively. These strongly suggest moderate to good 
porosity in the reservoir sands. The effective porosity 
values however reduced with noticeable differences with 
RSV 1, RSV 2, RSV 3, RSV 4 and RSV 5 recording values 
(measured in v/v) of 0.289, 0.252, 0.236, 0.236, 0.234 
and 0237 respectively. The differences observed between 
the density porosity and the effective porosity is 
interpreted to be microporosity contributions from 
shale/clay. Similar trend to this is also observed in Well B 
(Table III) and Well C (Table IV). 

2) Irreducible Water Saturation Swiry 
Observed irreducible water saturation is range from 0.09 
to 1 for Well A, 0.094 to 0.138 in Well B and 0.081 to 
0.087 in Well C (Table V). These values are high enough 
to ensure little to no water cut during production. 
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Table V Irreducible Water Saturation (Swirr ) in Langbodo 
Wells Reservoirs 


Reservoirs SwirrinWellA — Swirr in Well B Swirr in Well C 
RSV 1 0.088 0.094 0.081 
RSV 2 0.092 0.099 0.087 
RSV 3 0.100 0.101 0.081 
RSV 4 0.100 0.111 0.084 
RSV 5 0.092 0.138 0.083 


VII. CONCLUSION 


Careful comparison of the Tixier, Timur and Coates and 
Dumanoir equations and the mean (composite) permeability 
shows that the Timur’s relation closely agrees with the 
composite (mean) value, although some marked discrepancies 
were observed at some points. Petrophysically, subsurface 
reservoirs in the Langbodo field have reasonable hydrocarbon 
in their pore spaces, and estimated producibility indicators are 
good enough to support secondary migration of this oil into 
the borehole, if developed. Taking stand from the findings of 
this study, the use of a single method for estimating of 
reservoir permeability (K) is strongly discouraged. It is 
recommended that composite Permeability (K) equations be 
used. 
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